#!/usr/bin/env python

'''
program plots phase diagram for ferromagnetic triangular transverse
'''

import matplotlib.pyplot
import math
import numpy as np


#Last
B1 = [0,    0.05, 0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.7, 0.8,  0.9,  1.0, 1.1,  1.2,  1.3,  1.4,  1.5,  1.6,  1.7,  2.0]#, 3.0,   5.0, 6.0, 7.0, 10.0]
T1 = [2.21, 2.20, 2.20, 2.20, 2.18, 2.15, 2.12, 2.08, 2.04, 2.0, 1.95, 1.9, 1.85, 1.82, 1.81, 1.83, 1.86, 1.89, 1.94, 2.08]#, 2.72, 4.2, 5.1, 5.95, 8.35 ]

#middle
B3 = [0.1,   0.2,  0.3,  0.4,  0.5,  0.6, 0.7,  0.8,   0.9,   1.0,  1.1,  1.2,  1.3, 1.35  ]
T3 = [0.009, 0.02, 0.03, 0.04, 0.06, 0.07, 0.08, 0.105, 0.12, 0.15, 0.18, 0.21, 0.25, 0.28]

#smaller
B2 =        [0.05, 0.1,  0.2,  0.3,  0.4,  0.5, 0.6,  0.7,  0.8]
T2 =        [0.03, 0.09, 0.17, 0.25, 0.33, 0.42, 0.51, 0.62, 0.78]


matplotlib.pyplot.xlabel("B/J")
matplotlib.pyplot.ylabel("T/J")


matplotlib.pyplot.title('Position of the specific heat maximum. \n Triangular antiferromagnetic in transverse field. \nExact diagonalization')
matplotlib.pyplot.plot(B1, T1, "*-")
matplotlib.pyplot.plot(B2, T2, "*-")
matplotlib.pyplot.plot(B3, T3, "*-")

matplotlib.pyplot.show()